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We reexamine likelihood analyses of the Local Group (LG) acceleration, paying particular 
attention to nonlinear effects. Under the approximation that the joint distribution of the 
LG acceleration and velocity is Gaussian, two quantities describing nonlinear effects enter 
these analyses. The first one is the coherence function, i.e. the cross-correlation coefficient 
of the Fourier modes of gravity and velocity fields. The second one is the ratio of velocity 
power spectrum to gravity power spectrum. To date, in all analyses of the LG acceleration 
the second quantity was not accounted for. Extending our previous work, we study both the 
coherence function and the ratio of the power spectra. With the aid of numerical simulations 
we obtain expressions for the two as functions of wavevector and as- Adopting WMAP's 
best determination of as, we estimate the most likely value of the parameter j3 and its errors. 
As the observed values of the LG velocity and gravity, we adopt respectively a CMB-based 
estimate of the LG velocity, and Schmoldt et al.'s (f999) estimate of the LG acceleration from 
the PSCz catalog. We obtain (3 = 0.66±g;^; thus our errorbars are significantly smaller than 
those of Schmoldt et al. This is not surprising, because the coherence function they used 
greatly overestimates actual decoherence between nonlinear gravity and velocity. 

1 Introduction 

Comparisons between the CMB dipole and the Local Group (LG) gravitational acceleration 
can serve not only as a test for the kinematic origin of the former but also as a constraint on 
cosmological parameters. A commonly applied method of constraining the parameters by the 
LG velocity-gravity comparison is a maximum-likelihood analysis, elaborated by several authors 
(especially by Strauss et al. - -, hereafter S92). In a maximum-likelihood analysis, proper objects 
describing nonlinear effects are the coherence function (CF), i.e. the cross-correlation coefficient 
of the Fourier modes of the gravity and velocity fields, and the ratio of the power spectrum of 
velocity to the power spectrum of gravity. Here, with aid of numerical simulations we model the 
two quantities as functions of the wavevector and of cosmological parameters. We then combine 
these results with observational estimates of vlg and gLG> an d obtain the 'best' value of (3 and 
its errors. 

2 Modelling nonlinear effects 

We follow the evolution of the dark matter distribution using the pressureless hydrodynamic code 
CPPA (Cosmological Pressureless Parabolic Advection, see Kudlicki, Plewa & Rozyczka 1996, 
Kudlicki et al. 2000 for details). It employs an Eulerian scheme with third order accuracy in 



space and second order in time, which assures low numerical diffusion and an accurate treatment 
of high density contrasts. Standard applications of hydrodynamic codes involve a collisional 
fluid; however, we implemented a simple flux interchange procedure to mimic collisionless fluid 
behaviour. 

We studied the CF in an earlier paperP Here we use a different fitting function, which is 
more accurate for low values of k: 

C(k) = [l + (a k - a 2 k 15 + ai*; 2 ) 2 ' 5 ] ~°' 2 . (1) 

Parameters a/ were obtained for 35 different values of as in the range [0.1,1], and we found the 
following, power-law, scaling relations: 
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The fit was calculated for k E [0, 1] h/Mpc, with the imposed constraint C(k = 0) = 1. 

We have found that the ratio of the power spectra, 1Z, obtained from simulation can be fitted 
with the following formula: 

1l(k) = [l + (7.071fc) 4 ]" Q , (3) 

with 

a = -0.06574 + 0.29195cr 8 for 0.3 < cr 8 < 1 . (4) 
3 Parameter estimation 

We apply our formalism to the PSCz survey. As the value of as we adopt its WMAP's estimate, 
as = 0.84 (±0.04; Spergel et al. 2003). This specifies the coherence function and the ratio of the 
power spectrum of velocity to the power spectrum of gravity. Also, this provides a normalization 
for the power spectrum of density. 

3.1 Parameter dependence of the model 

The likelihood of specific values of (5 and of the linear bias, b, is determined by the following 
distribution (Juszkiewicz et al. 1990, Lahav, Kaiser & Hoffman 1990): 
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where a s and a v are respectively the r.m.s. values of a single Cartesian component of gravity 
and velocity, (x, y) = (g/fj g , v/cr v ), and tp is the misalignment angle between g and v. Finally, r 
is the cross-correlation coefficient of g m with v m , where g m (v m ) denotes an arbitrary Cartesian 
component of g (v). 

In this distribution, the observables are g and v, or g, v, and the misalignment angle, 
-0. Following Schmoldt et al.^ (hereafter S99), we adopt for them the following values: g = 
933 km • s _1 (from the distribution of the PSCz galaxies up to 150 h~ l Mpc), v = 627 km • s _1 
(inferred from the 4- year COBE data by Lineweaver et al. 1996), and tjj = 15°. 

The theoretical quantities are a s , a v , and r. The variance of a single spatial component of 
measured gravity, a 2 ,, is a sum of the cosmological component, a^ c , and errors, e 2 . Since gravity 
here is inferred from a galaxian, rather than mass, density field, we have <t 2 c = b 2 Sg, where s 2 
is the variance of a single component of the true (i.e., mass-induced) gravity, 

i r°° - 
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Figure 1: Likelihood contours for the parameters (3 and b, corresponding to the confidence levels of 68, 90 and 
95%. The maximum of the likelihood function is denoted by a dot. 



The gravity errors are twofold: due to finite sampling of the galaxy density field, and due 
to the reconstruction of the galaxy density field in real space. Therefore, e 2 = (<t| n + <7 2 ec )/3, 
where Og N and respectively the shot noise (or, sampling) variance and the reconstruction 

variance. In brief, 
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Finally, errors in the estimate of the LG gravity do not affect the cross-correlation between the 
LG gravity and velocity, but increase the gravity variance. This has the effect of lowering the 
value of the cross-correlation coefficient. Specifically, we have 
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Thus, the likelihood depends explicitly on the parameters b and J7 m , or on b and (3 = f2^ 6 /&. 



3. 2 Joint likelihood for [3 and b 

Figure ^ shows isocontours of the joint likelihood for j3 and b, corresponding to the confidence 
levels of 68, 90 and 95%. The maximum of the likelihood is denoted by a dot. The corresponding 
values of (3 and b are respectively 0.66 and 0.94. The isocontours of the likelihood are much 
more elongated along the 6-axis than along the /3-axis, making the resulting constraints on £l m 
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Figure 2: The marginal distribution for j3. The result of marginalizing over all possible values of b (from zero to 
infinity) is shown as a dashed line. The result of marginalizing over the values of b in the range [0.7, 1.1] is shown 

as a solid line. 

much weaker than on [5. In practice, therefore, from our analysis we cannot say much about 
bias, or TO , alone. On the other hand, we can put fairly tight constraints on (3. 

To do this, we use the fact that the square root of the variance of the PSCz galaxy counts 
at 8 h- 1 Mpc is a^ SCz ~ 0.75 (Sutherland et al. 1999). Combined with WMAP's estimate of 
<Tg, this yields for the bias of the PSCz galaxies the value about 0.9. Therefore, we adopt here a 
conservative prior for the bias, namely that it is constrained to lie in the range [0.7, 1.1]. These 
limits are marked in Figure ^ as dashed vertical lines. We marginalize the likelihood over the 
values of b in this range. The resulting distribution for j3 is shown in Figure |2] as a solid line. 
We obtain (3 = 0.66±8$ (68% confidence limits). 

4 Summary and conclusions 

The analysis of the LG acceleration performed by S92 and S99 was in a sense more sophisticated 
than ours. Both teams analysed a differential growth of the gravity dipole in subsequent shells 
around the LG. Instead, here we used just one measurement of the total (integrated) gravity 
within a radius of 150 h^ 1 Mpc. Nevertheless, the errors on (5 we have obtained are significantly 
smaller than those of S92 and S99. In particular, S99 obtained j3 = 0.70 j^;^ at la confidence 
level. It is striking that while our best value of (3 is close to theirs, our errors are significantly 
smaller. The reason is our careful modelling of nonlinear effects. In a previous paper we showed 
that the coherence function used by S99 greatly overestimates actual decoherence between non- 
linear gravity and velocity. Tighter correlation between the LG gravity and velocity should 
result in a smaller random error of f3; in the present work we have shown this to be indeed the 
case. 

Acknowledgments 

This research has been supported in part by the Polish State Committee for Scientific Research 
grants No. 2.P03D.014.19 and 2.P03D.017.19. 

References 



1. M.A. Strauss, A. Yahil, M. Davis, J.P Huchra and K. Fisher, ApJ 397, 395 (1992) (S92) 

2. M.J. Chodorowski and P. Ciecielag, MNRAS 331, 133 (2002) 

3. I. Schmoldt, et al., MNRAS 304, 893 (1999) (S99) 



